July 19, 2010 0:43 WSPC/INSTRUCTION FILE paper 



O 

(N 



o 
o 



International Journal of Modern Physics D 
© World Scientific Publishing Company 



Ellipticity of Structures in CMB Sky Maps 



Ralf Aurich, Holger S. Janzer, Sven Lustig 

Institut fiir Theoretische Physik, Universitdt Ulm, 

Albert- Einstein- Allee 11, D-89069 Ulm, Germany 



o 

^^ ' Frank Steiner 

,il^ ■ Universite Lyon 1, Centre de Recherche Astrophysique de Lyon 

Qh' CNRS UMR 5574, 9 avenue Charles Andre, F-69230 Saint-Genis- Laval, France 

6 

-^-J , We study the ellipticity of contour lines in the sky maps of the cosmic microwave back- 

~3 ' ground (CMB) as well as other measures of elongation. The sensitivity of the elongation 

I I ' on the resolution of the CMB maps which depends on the pixelization and the beam 

profile of the detector, is investigated. It is shown that the current experimental accuracy 
does not allow to discriminate between cosmological models which differ in curvature by 
^ ' A^tot = 0.05. Analytical expressions are given for the case that the statistical properties 

fvj ' of the CMB are those of two-dimensional Gaussian random fields. 
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f--^ . 1. Introduction 
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The cosmic microwave background (CMB) gives a wealth of information about the 
Universe. The usual statistical analysis comprises the angular auto-correlation func- 
tion C{d) as well as the angular power spectrum 5Tf ^ which lead to an improved de- 
termination of the cosmological parameters. Besides these statistics there are topo- 
logical/geometrical descriptors of the CMB such as the number densities of maxima 



C^ ' and minima, the ellipticities of the peaks, and peak correlation properties. Under 

the assumption that the CMB has the statistical properties of two-dimensional ho- 
mogeneous and isotropic Gaussian random fields, the theoretical predictions are 
derived in Ref ill 

This paper deals with the structure of the contour lines of a given temperature 
level 5T in the CMB maps. In the vicinity of extrema these contour lines can be 
approximated by ellipses having the ellipticity E 

E := ^ , E > 1 , (1) 

where a > b denote the major and minor semi-axes of the ellipses which best match 
the contour line around a given peak. The interesting question is how the ellipticities 
are distributed for a given map in dependence on the temperature level ST. This 
provides a further tesl!^'^ of Gaussianity by comparing the ellipticities with the cor- 
responding results of two-dimensional Gaussian random fields '. In Refs.|4l[5l|6]it is 
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shown, however, that there are better geometrical estimators as tests of Gaussianity 
than the ellipticity like the Gaussian curvature. Therefore, the main motivation for 
this paper is to study the discrepancy of the ellipticity between the observations 
and simulations of the CMB as reported in Refs. [H [3 [9l [TOl [HI [HI In these papers 
an excess in the ellipticity E = 2.2 . . . 2.7 is found compared to standard ACDM 
models having E around 1.65. Due to noise and the foreground uncertainties it is 
difficult to obtain the ellipticity from the measured CMB radiation. The ellipt icity 
was determined from the COBE-DMR ma^ and from the BOO MERa nG mapPSI. 
These results are confirmed by the analysis of the WMAP lyr dat a I ^and the 3yr 
datsM^. Furthermore, a strong dependence of the ellipticity on the curvature of the 
Universe is claimed and, with respect to the excess ellipticity, this is interpreted as 
a hint towards a hyperbolic Universe. In Ref.|4] there is no such ellipticity-curvature 
correlation found by studying low fl and flat universes based on CDM models. Thus 
there remains the question whether the excess ellipticity is real and how such an 
excess, if present, has to be interpreted. 

2. CMB Sky Maps 

2.1. Cosmological models 

The standard AGDM model of cosmology is given by a certain set of cosmological 
parameters. The curvature of the universe is determined by the value of the total 
energy density fitot at the present epoch which is the sum 

f^tot = ^A + f^cdm + r^b + ^r (2) 

of the energy density f^A of dark energy, ilcdm of cold dark matter, fib of baryonic 
matter and il^ of radiation. A value fitot > 1, = 1 or < 1 reveals a spher ical, flat 
or hyperbolic universe, respectively. The concordance ACDM mode l I I leads to 
a flat universe where the angular power spectrum 6T^ — l{l + l)C;/(27r) has the 
flrst acoustic peak at I ~ 220. The following analysis investigates the dependence 
of the ellipticity on the curvature of the universe. To that aim only the amount of 
dark energy J7a is varied, all other cosmological parameters are those of the ACDM 
concordance model. 

In additio n to th e infinite volume ACDM concordance model, a model with a 
cubic topolog y I I is also studied which has a finite volume and is statistically 
anisotropic. This model thus violates some conditions required for the following an- 
alytical expressions as outlined in the Appendix. This model serves as a test whether 
an anisotropic non-trivial topology can be discerned by an elongation measure. The 
sky maps are simulated for the cubic topology using the same cosmological parame- 
ters as for the infinite concordance model. The side length L of this toroidal topology 
is chosen as L = 3.86 Lh where Lh is the Hubble length. For this length L a better 
agreement with the correlations observed in the CMB sky is found '-' ^" ^ ' than for 
the ACDM concordance model. The eigenmodes belonging to the first 50 000 eigen- 
values for the cubic topology are used for the simulation, i. e. a total of 61 556 892 
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eigenmodes. Furthermore, the eigenmodes are expanded in the spherical basis up 
to Imax = 1000 yielding sky maps with structures beyond the third acoustic peak 
where the Silk damping already smoothes a further fine-structure. In this way CMB 
sky simulations are obtained for which resolution effects can be studied. 

For a given cosmological model, i. e. a given set of cosmological parameters, one 
can calculate the multipole spectrum 

c^--=(^,tM') > (3) 

\ m— — l I 

where aim are complex coefficients obtained from the expansion of a CMB sky map 
5T{6,<p) into spherical harmonics Yi^{<p^6) due to 

oo I 

6T{e,cl>) =Y.Y1 ^lmYl,n{<P,9) . (4) 

The multipole spectrum Ci is obtained from an ensemble average denoted by (. . .) 
over infinitely many realisations of universes with fixed cosmological parameters. 
The deviation of the multipole spectrum C; of an individual realisation from the 
ensemble average is characterised by the cosmic variance 

Var(G) := ^ , (5) 

where one has to assume that the CMB is a homogeneous isotropic Gaussian random 
field. This has to be taken into account by comparing theory and experiment. 

2.2. Resolution of the sky maps 

Statistical measures of the niveau lines depend sensitively on the resolution of the 
sky maps. Thus a few remarks are in order. A natural cut-off in the multipole 
space is provided by the physics of the CMB, especially by the Silk damping and 
the smoothing due to the thickness of the surface of last scattering. One would 
need sky maps with a resolution of at least this physical cut-off in order to capture 
all genuine CMB structures. This is currently beyond the possibilities, and one is 
forced to consider sky maps which are limited by the measurements, i. e. by the 
beam profile of the detector. The comparison of the measured sky map with a 
simulated map requires that the simulation is accordingly smoothed. A symmetric 
Gaussian smoothing kernel is sufficient for most applications which is a special case 
of a general symmetric smoothing kernel represented by the window function Fi. 
The smoothing operation is done in multipole space by 

oo I 

5T{e,4>) ^ 5Tf,{0,^) =Y.Y. FiaimYi^{e,c^) . (6) 

(=0 m=-l 

The symmetric Gaussian kernel is given by the window function 



^Gauss ^ g^p 



(7) 
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where cTg is the width of the symmetric Gaussian kernel which is usually parame- 
terised by its full width at half maximum tTfwhm- The conversion formula is given 

by 

TT f^fwhm 



(8) 



* 180° 2V2M2) ' 

where Ufwhrn is given in degrees. 

An analysis of the temperature fluctuation field in position space requires a 
pixelization of the data on the sphere. The CMB sky maps are usually discretized 



in the HEALPix^^ format where every pixel covers an equal area. The resolution 
parameter A'^gidc defines the total number of pixels by N^^ = 12 N'^^^^- The resolution 
in position space can be defined by the square root of the area of a single pixel. The 
area of a pixel is 10800/ (ttA^^.^ J deg^ 

The WMAP data are available in N^ide = 512 at the LAMBDA website 
(lambda.gsfc.nasa.gov), and the Planck data will be provided in A^side = 2048. As a 
rule of thumb the pixel resolution should relate to the resolution in multipole space 
by In 



2.3. Contamination of CMB measurements 

A measurement gives the genuine CMB signal superimposed with the emission of 
foreground sources and the detector noise. 

Noise is unavoidable in measurements of CMB sky maps and acts mostly on the 
smallest scales, i. e. at large multipoles I. It can be modulated by generating in each 
pixel random fiuctuations which are added to the simulated CMB signal. In the case 
of the WMAP data the standard deviation of the random fiuctuations is propor- 
tional to l/VA'obs where the constant of proportionality is stated on the LAMBDA 
website and A^obs is the number of observations of a given pixeJ^. Since A^obs de- 
pends on the pixel, the noise also depends on the direction yielding anisotropic noise 
properties. 

Astrophysical radiation sources bring foreground contaminations on measured 
CMB sky maps. Cleaning operations reduce their contributions, but cannot avoid 
that still some sky regions have to be excluded from the analysis by masking them 
out. This is a serious problem for estimators which are based on full sky informa- 
tion as the angular power spectrum ST^. For the structure analysis the remaining 
foreground contamination in regions, that are not masked out, is more important 
since it changes the properties of the niveau lines. 



3. Measures of Elongation 

Probably the most common measure of the elongation of a structure is the ellipticity 
E which is based on the assumption that contour lines are approximated by ellipses. 
E is determined by the semi-axes a and 6, which are real numbers with a > b, oi the 
best-fit ellipsis. To parameterise the elongation, various dimensionless combinations 



July 19, 2010 0:43 WSPC/INSTRUCTION FILE paper 



EUipticity of Structures in CMB Sky Maps 5 



of the semi-axes can be used. Three quantities are found in the hterature and given 
by the elhpticity E, the eccentricity e, and a further efongation measure e. Their 
definitions and conversion formulas are given by 



E 



l + 2e 
1 -2e 



1 



VT^l^ 



l< E < oo 



(9) 



V^^^¥ I 1 „ 



l + 2e 



0< £ < 1 



and 



(10) 



1 1 



2 1 



b2 



2 2-e^ 



1 1 



1 



2 1 



1 



0< e < 



(11) 



Larger values indicate a larger degree of elongation for all three measures. There are, 
however, further possibilities to quantify structures, see Section[5]and the Appendix. 



4. Cosmological Dependence of Gaussian Random CMB Maps 

4.1. Number of maxima and minima 

Before we turn to the elongation at maxima and minima, we will discuss the number 
-^max and Nmin of maxima and minima per solid angle on a 2-sphere. Assuming 
that the temperature fluctuations 6T (fi) are given by a homogeneous and isotropic 
Gaussian random field, the corresponding derivations of the analytical formulae for 
the ensemble average of the number of maxima and minima per solid angle on a 
2-sphere are given to some extent in Ref. |T, and further relations can be found in 
[Appendix A.l 

In the case of a homogeneous and isotropic universe the distribution of -/Vmax as 
a function of the temperature threshold ST is determined completely by the three 
parameters CTq, af and cr|, which are given by 

(/ + n)! 



E 



2/ + 1 

An 



Ci\Fif 



{l-n)\ 



with 



n = 0,1,2 



(12) 



as discussed in the Appendix (see (|A.2|) . (jA.5|) and (|A.6|) ). Obviously these param- 
eters can be calculated from the power spectrum of the CMB. In this subsection 
we consider the number iVmax of maxima per solid angle depending on the nor- 
malised temperature v := — where cto is the standard deviation of the temper- 
ature fluctuations. The corresponding number of minima is given by the relation 

Nmin{v) = iVmax(-i^)- 

For a homogeneous and isotropic Gaussian random field, the distribution 
Nm?,x{v) is given by 



A^max(i^) 



1 



(27r)t 



5*2 



exp 



G(z/,7,a) 



(13) 
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Fig. 1. The ensemble average of A^maxCi^) for the ACDM concordance model computed by Eq. II13I I 
is shown as a solid line. Panel (a) compares this curve with the average of Afmax(i') computed 
from 50 realisations of the infinite volume model (circles). Panel (b) shows the result obtained 
from 50 realisations of the torus universe with the side length L = 3.86 Ljj (circles). To reveal the 
cosmic variance, the la standard deviation is displayed as a grey band which is calculated from 
the corresponding 50 maps. 



with 



G(i/,7,a) := 71^(1-7^) 



exp 



'2(1-7^) 



V2vr(l-72) 



2 W2(l-7^), 



exp 



l+2a^(l-7^) 



^2a2(l-72) + l 



1 _ ierfc I — 

2''''l V2(l-72)(l + 2a2(l-72)) 



and is derived in [Appendix A.l see Eq. (|A.23|) . Here erfc(a;) is the complementary 
error function. The parameters are 



a 




1 



2a? 
'2 



and 



CTo 



7 



(72 (To a 



(14) 



In Fig. [T] the theoretical distribution (jlSp is shown and compared with the mean 
value obtained from 50 CMB maps simulated at the HEALPix resolution A'sidc = 
512 using a smoothing of 1° and I < 1000. Both panels are based on the cosmological 
parameters of the infinite volume best-fit ACDM model of the WMAP data. In 
panel (a) the infinite model, which is isotropic, is shown, whereas panel (b) shows a 
multi-connected model which is anisotropic and thus does not fulfil the assumptions 
on which (fT3|) is based. However, as revealed by Fig. [ijb) the data points match 
the distribution (fT3|) with the same quality as in Fig. [TJa). In addition, the 50 
realisations can be used to get an estimate of the cosmic variance. The corresponding 
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la standard deviation is displayed as a grey band in Fig. [TJ The results for the 
mean value and the variance of iVniax('^) for the torus universe are almost the 
same as for the infinite one. One could have expected this result since the number 
of maxima is primarily determined by the smallest scales, and a multi-connected 
model, which possesses in the shown case a cubic topology with a side length L — 
^.86 Lh, differs from the infinite one only on very large scales. Similar results are 
also expected in the case of other topologies, homogeneous and inhomogeneous 
manifolds, if the topological scales of the universe are of comparable order. Only 
if the volume of the universe is much smaller than the volume inside the surface 
of last scattering, one can expect any differences in the iVmax distribution to those 
calculated in |Appcndix A.l| In such a case one would get an influence of the topology 
on small scales. But until now no hint is found in the data for a model with such a 
sufficiently small volume. 
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Fig. 2. In panel (a) the formula of the total number A^max of maxima per solid angle resulting 
from Equation (A1.7) in Ref. [T] (solid line) is compared to our Eq. 11511 (dashed line). In panel 
(b) the parameter a is plotted as a function of the smoothing scale o-f^hm '^°'' three cosmological 
models. All values of a are very close to one. 



Our formulae for the various A'max distributions in [Appendix A. l] differ by terms 
depending on the parameter a from the corresponding equations in Ref. [T] Here 
these differences are discussed using the example of the total number of maxima 
per solid angle. Our formula to this quantity is given by 



iVn 



= N, 



B&E 



[(g^ - 1) VI + 2q2 + 1] a/3 
Vl + 2a2 



(15) 



which is derived in [Appendix A.l see (JA.19P . The corresponding result obtained 
from Equation (Al.7) in Ref. 1 is N^^^ = ^^J"^ y- , which is the leading term of the 



Laurent series of our Eq. (fT5|) at a^ = 1. Both formulae are compared in Fig. [2Ja) 



where iirO* iVmax is displayed. In the case of A^max ^^ a- independent, i.e. model 
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independent value of -y= is obtained. In contrast our result for AttO*'^ -/Vmax depends 
on a and therefore on the cosmological model since a, al and cr| depend on the an- 
gular power spectrum. In Fig.[2Ja) we have varied a € [1, 2] in order to demonstrate 
the difference between iVmax and the leading term N^^^. Is this interval of a re- 
alistic for cosmological models compatible with the measured data? To answer this 
question we have calculated af and cr| with a cut-off at I — 1000 with the angular 
power spectrum of the best-fit ACDM model of the WMAP data. In addition, these 
parameters are computed from the angular power spectrum of two models with 
positive and negative curvature, where the cosmological constant is varied and the 
other cosmological parameters are held fixed. These three cosmological models lead 
to different values of a, whose dependence on the smoothing scale is displayed in 
Fig.mjb). One observes that a is almost equal to 1, and thus 9*^ is very small for all 
realistic cosmological parameters and smoothing scales. Therefore, we conclude that 
for all practical purposes the additional dependence of 4it9*'^ -/Vmax on a is without 

distributions in Ref. |T| are adequate. 



4.2. Elongation at Maxima and Minima 

Now we turn to the elongation measure e at maxima and minima. The ensemble 
average and the corresponding second moment of this quantity are given by 



2a^ + 3-|^/^^ln 



1+^ 



-^ 



4[(a2-l)Vl + 2a2 + l] 



(16) 



and 



1 - 3 VI + 2a2 + 4(V1 + 2a2 _ i) (^^) 

which both depend on the parameter a. The formulae are derived in the 
[Appendix A. 2 The mean value of the elongation e together with the standard devi- 



ation due to the cosmic variance is displayed in panel (a) of Fig. [S] The distribution 



Pie) = 24e(l-4e2)a^v/rT2^ ^^^^ 

(1 + 8a2e2)* [(a2 - l)Vl + 2^2 _^ ^ 

of the elongation is plotted for four values of a in Fig.[3l^b). This demonstrates that 
P(e) depends on the parameter a and thus in turn via the angular power spectrum 
on the cosmological parameters. In this way the elongation measure e encodes some 
properties of the Universe. 

This dependence of the elongation e on cosmological parameters is not discussed 
in Ref. [1] In currently admissible cosmological models the dependence on a in 
(jA.30|) can be neglected since the corresponding values of a are very close to one. 
The ensemble average of the elongation e is approximatively given by 0.197 for 
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Fig. 3. In panel (a) the ensemble average of the elongation e with the corresponding la standard 
deviation is displayed depending on the parameter a. In panel (b) the distribution P{e) of the 
elongation e is plotted depending on the parameter a. 

a = 1. Similar results are obtained for the other elongation measures, i.e. for the 
ellipticity E and the eccentricity e. The ensemble averages of these quantities are 
given by (|A.36p and (|A.33|) . Their values are approximately E « 1.648 and e « 0.715 
independent from the model. The cosmological influence on the elongation e, the 
ellipticity E or the eccentricity e is very difficult to measure. The difficulty arises 
from the large cosmic variance for these quantities which is displayed as a grey band 
in Fig. ^a.) in the case of the elongation e. Furthermore, also the detector noise 
and the foregrounds in a sky map obtained by observations have an influence on 
the elongation, but this will be discussed later. 

In addition, it is not expected that the topology of the Universe, i. e. a multi- 
connected spatial space, can influence the elongation of contour lines since the 
volumes of the fundamental cells are too large. Nevertheless, several measures of 
the elongation are discussed below for the special case of the cubic topology which 
is motivated by the fact that this model is statistically anisotropic. 



5. Elongation of Hot and Cold Spots 

There are different algorithms to compute the ellipticity. One method is based on 
the Taylor expansion of the temperature field at local maxima or minima as used in 
the Appendix. Alternatively, the ellipticity can be calculated by using an "inertia" 
tensor as described in the following. 

5.1. Computation of the ellipticity 

In the following the normalised temperature field 

STp, (h) 
u{n) = '^ ' (19) 
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is used. In the theoretical derivation which is outlined in the Appendix, this ctq is 
obtained by the ensemble average by using uq — ^/CJfi) where C(0) is the 2-point 
correlation function of the CMB at a separation angle d — Q. Such an ensemble 
average is not possible for a CMB sky map obtained by observations, of course, 
since there is only a single CMB sky from our point of view. Therefore, a sky map 
analysis computes the value of ctq from a given single sky realisation as the usual 
standard deviation of the temperature field 



^^^\Nsn / {6TF,{n)~fYdn , (20) 

y Jsn 

where T is the mean temperature. In the analysis of masked maps the observed 
area is (557 and Nsq is the fraction of the sky that is not masked. 

To get a binary image we define the excursion set of hot spots (HS) by 

qHS ^ {neS^\u{h) >i^} . (21) 

Here one excludes all regions with a temperature value lower than the value of the 
threshold ly. Analogously the excursion set of cold spots (CS) is defined by 

Q[r^ - jneSXn) <u} . (22) 

The multi-connected set Q^, can be interpreted as the union 

Q. = [JQl (23) 

i=l 

of the Ni, simply connected spots Q]^. We compute the ellipticity for every spot Q\, 
by using the inertia tensor defined by 

rCOM — ( -^^^ Ixy \ _ ( JJy^dxdy ~J Jyxdxdy \ 

^^''' ■" I ly. lyy J \ - f f xy dxdy JJx^dxdy J • ^^^^ 

The spots are projected onto the xy-plane where the origin of the coordinate system 
matches the centre of mass of the spot Q^. The symmetric tensor I{Ql,) has two 
real eigenvalues Ai and A2 whose relation to the ellipticity E is given by 

E{Q\^ = \[^ . (25) 

V M 

For the special case of an ellipse with axes a and b the eigenvalues can be computed 
to be 

Al = -a^b and A2 = -ab^ with Ai > A2 . (26) 

4 4 

Inserting this into (PSJ) leads back to the definition (HI), respectively ^ of the 
ellipticity 

V A2 
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The other quantities for elongation given in (|10p and (ITlT) can be obtained from the 
conversion formulae. The true niveau lines are not perfect ellipses, of course, but 
the above prescription can be used in order to obtain the ellipticity of that ellipse 
which approximates best the niveau line. Since there are niveau lines which deviate 
strongly from an ellipse one needs a further criterion to eliminate such curves in 
order to avoid non-sense results. 



5.2. The dependence of the ellipticity on the resolution 

The ellipticity depends sensitively on the accuracy of the sky map and care has to 
be taken with respect to the resolution of the map, the beam profile and the noise 
properties. In so far, it is not a very robust measure to characterise the CMB sky, 
and the details of the computation of the ellipticity have to be defined clearly in 
order to obtain reproducible values. 

Let us at first address the issue of the resolution of the sky map. The CMB sky 

1 S 
maps are usually stored in the HEALPix format whose resolution is determined 

by the parameter Nsidc- In order to determine the dependence of the ellipticity 

E on the value of A^side, a sky map is simulated for the cubic topology using the 

cosmological parameters of the ACDM concordance model. We start with a highest 

resolution of A'sidc — 4096. The sky maps are normalised such that the temperature 

fluctuations possess a standard deviation of 66/zK, which is the value obtained from 

the five-year ILC map outside the KQ75 mask. 

For such a simulation the ellipticity E is computed for the contour lines of a 
given temperature threshold 5T, and the mean value is shown in Fig. |4l Contour 
lines that are too small are excluded since they are determined by very few pixels 
such that the ellipticity E is ill defined. In Fig.|4lja) the criterion is that the contour 
lines should at least circumference 10 pixels with respect to the chosen resolution 
iVsidc. This implies that ever more small contour lines are taken into account as the 
value of Aside increases. Thus the mean values are computed from different sets of 
contour lines. This contrasts to the criterion of a fixed minimal area which leads 
to the result displayed in Fig. |4ljb). There, all contour lines are taken into account 
which are larger than 10 pixels with respect to the A'sidc — 256 resolution, i. e. all 
contours encompassing an area larger than 0.525 deg . Thus the curve belonging to 
Aside — 256 is the same in both panels. Whereas the exclusion criterion of a fixed 
pixel number leads to a robust estimation of E only at very large values of A"sidc, the 
exclusion criterion of a fixed area gives consistent values already for values of A'side 
as small as 512. But recall from Section 221 that the expected ellipticity is i? ~ 1.648 
for a Gaussian field. This value is, however, obtained with the first selection criterion 
which takes an increasing number of contour lines into account. This analysis thus 
favours the first selection criterion and, furthermore, demonstrates the sensitivity 
of the elongation measures. 

The dependence of the ellipticity on the smoothing is not very strong. This 
can be inferred from Fig. El^b) where the value of a is plotted as a function of the 
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Fig. 4. The ellipticity i? is computed for a realisation of the cubic universe as described in the 
text based on the cosmological parameters of the concordance model. The dependence on the value 
of A'gjjg is presented. Panel (a) shows the ellipticity E where all contours are taken into account 
that enclose at least 10 pixels. The ellipticity E increases with decreasing N^^^^. In Panel (b) all 
contours are selected enclosing an area of at least 0.525 deg. 



smoothing crfwhm- Up to a smoothing of 1 degree all values of a are below 1.00004 
and the expectation value for the ellipticity E, Eq. (|A.36p . is nearly constant for 
these smoothings. This behaviour is confirmed in Fig.[5]where different smoothings 
are applied to the sky maps with a fixed HEALPix resolution of A'sidc — 4096. 
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Fig. 5. For the CMB simulation used in Fig.|4{a) having the HEALPix resolution of Afgj^p : 
the ellipticity E is shown in dependence on the width of the Gaussian beam. 
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5.3. The dependence of the elongation on the curvature 

Let us now turn to the dependence of the elongation on the curvature which is 
claime d I I I I I I to be sufficiently strong in order to reveal the curvature of 
the Universe. However, in Ref. 3] no such ellipticity-curvature correlation is found. 
Since the elongations depend on the parameter a and thus on the cosmological 
parameters, such a correlation could exist. The discussion in the previous section 
has shown that the elongation depends on the size of the niveau lines which have 
to enclose at least n^nin pixels. For the computation of the elongation our algorithm 
requires at least an area greater than 4 pixels. The number of niveau lines which 
can be used in a statistic depends on the normalised temperature v. It is maximal 
at z^ = ±1 where the number N^, of spots has its maximum. The values of n,nin 
and of the normalised temperature v are required to specify the mean ellipticity 

The ensemble average of the elongation and the corresponding cosmic variance 
are computed from an ensemble of 1000 sky realisations. We analyse the excursion 
sets Qi, for hot spots with thresholds i' > and for cold spots with thresholds i' < 0. 
We get for every threshold v an amount of N^, spots with individual areas f^((^^) 
and elongations e(Q1). The mean ellipticity E{v,naihi) is computed at i/ = 1 by 
averaging the values of all spots with n > n,„in. This method is used in Refs. [Qj [TOl 
111) 1121 An alternative averaging is used for the mean elongation e(i^, n, An) which 
is computed as a moving average around a spot size n, i. e. from spots having pixel 
numbers within the interval [n — An,n + An]. Our examination shows that these 
mean values are very stable with respect to the choice of the threshold i^ in a range 
of i/e [-2,2]. 
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Fig. 6. The ensemble average of the mean eUipticity E(v, nuiin) with v = 1 (left) and the mean 
elongation e{y, n, An) with u = 1 and An = 20 pixels (right) is computed using three Gaussian 
smoothing kernels (from top to bottom: Cf^jjuj = 20 arcmin, 40 arcmin, and 60 arcmin). The grey 
bands reflect the cosmic variance of the flat ACDM concordance model. The ensemble average is 
computed from 1000 realisations. The sky maps have a resolution of Afgjjg = 512. 
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In Fig.|6]the ensemble average of the mean ellipticity E{i', n„^in) with v — 1 (left) 
and the mean elongation e{i>, n, An) with v — 1 and An — 20 pixels (right) is shown. 
The average is computed from 1000 realisations. The sky maps possess a resolution 
of iVgide = 512 and /max — 1000. The grey band reflects the cosmic variance of 
the flat ACDM concordance model (solid line). This provides the criterion whether 
two cosmological models can be distinguished by an elongation measure. In order 
to distinguish them the corresponding curves should deviate more than the width 
of the band due to the cosmic variance. Note that in contrast to the previous 
figures, the Fig. [5] displays the elongation as a function of the spot size and not 
of the temperature to which the niveau lines belong. This kind of plotting reveals 
the remarkable fact that larger spots possess a higher degree of elongation than 
smaller ones. This dependence is not covered by our analytical formulae which 
are based on the infinitesimal neighbourhood at maxima or minima. The impact 
of three Gaussian smoothing kernels is also investigated in Fig. |6] (from top to 
bottom: Cfwiim — 20,40,60 arcmin). An increased smoothing reduces the degree 
of elongation. This is due to the symmetric smoothing kernel which rounds off the 
structures. Here the resolution of A^'sido = 512 does not resolve all physical structures 
of the Gaussian temperature field. This contrasts to Fig. \5\ which is based on a sky 
map with a resolution of A'sido = 4096 being sufficiently fine grained to resolve all 
structures such that the analytical formulae are applicable. 

In Fig.[n]the ensemble average of the elongations is plotted for three cosmological 
models which differ with respect to their curvature. In addition to the fiat ACDM 
concordance model, the ensemble average of a positively curved universe (dotted 
line) and a negatively curved universe (dashed line) is displayed. The figure reveals 
that a universe with negative curvature possesses a larger degree of elongation for 
a given spot size than the fiat universe. This trend is continued for a universe with 
positive curvature which possesses a smaller degree of elongation for a given spot 
size than the flat universe. This encouraging behaviour could be used to distinguish 
different cosmological models if their corresponding curves differ by more than the 
scattering due to the cosmic variance. It turns out that this depends on the reso- 
lution of the sky maps. By applying wider smoothing kernels (from top to bottom 
in Fig. |6]) the differences between the curves decrease and are getting insignificant 
compared to the cosmic variance for smoothings around crfwhm = 60 arcmin. Note 
that the mean elongation e{i>, n, An) based on a moving average does a better job 
than the mean ellipticity E{iy,nyain)- This is caused by the fact that the differences 
in the elongation vanish for spots with a large area, and thereby their inclusion 
blurs the signal. The figure reveals that a resolution of at least crfwhm = 20 arcmin 
is necessary in order to distinguish between cosmological models having a difference 
in curvature of Ailtot — 0.05. This excludes the application of the ILC sky map 
of the WMAP team which has a resolution of about Cfwhrn = 60 arcmin. Better 
resolved sky maps such as the W-band maps of the WMAP team are available but 
they are contaminated with noise. Hence we now focus on the influence of noise on 
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Fig. 7. The ensemble average of the mean eUipticity Eiy, n^j^) with i^ = 1 (left) and of the 
mean elongation e(!/, n, An) with i/ = 1 and An = 20 pixels (right) is presented. These curves are 
obtained from 1000 sky simulations which take the window function of the W-band (channel 4) of 
WMAP into account. The grey bands reflect the cosmic variance of the flat ACDM concordance 
model (solid curve). The three rows give the results without noise (top), with noise (middle) and 
a comparison of both cases for the flat ACDM concordance model (bottom). The sky maps have 
a resolution of Af^j^g = 512. 
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In order to address the question whether the W-band maps are suitable for 
an elongation analysis, 1000 sky maps are simulated which tak e in to account the 
explicit window function of the W-band (channel 4) of WMAP^. In Fig. [7] the 
ensemble averages obtained from such sky maps are shown for the mean ellipticity 
Eiy, nniin) with v = 1 (left) and for the mean elongation e(;^, n. An) with v = 1 
and An = 20 pixels (right). Again the grey band reflects the cosmic variance of the 
flat ACDM concordance model (solid line) . In addition, a positively curved universe 
(dotted line) and a negatively curved universe (dashed line) is displayed. 

The upper row of Fig. [7] neglects noise completely and only takes the W-band 
window function into account. Both figures are similar to the two upper figures 
in Fig. [S] which have a resolution of (Tfwhm — 20 arcmin comparable to that of 
the W-band. The mean elongation e{v^ n, An) of spots with a small area reveals a 
saturation which is absent in Fig. [6l Furthermore, both elongation measures yield 
smaller values for spots with a larger area. This is caused by the application of the 
W-band window function which deviates strongly from a Gaussian. In the middle 
row of Fig.[7]the noise is taken into account which changes the behaviour drastically. 
No useful discrimination between the different curved models is possible any more. 
The lowest row in Fig. [7] shows a comparison for the ACDM concordance model 
with and without noise in one figure in order to emphasise the deteriorating effect 
of noise. A further analysis shows that an incomplete sky coverage additionally 
broadens the cosmic variance. Therefore, no comparison to measured data can be 
shown here. 

Finally, we compare the infinite volume ACDM model with a statistically 
anisotropic multi-connected model for which we again choose the cubic topology 
with side length L = 3.86 Lh- As already discussed in Section H?T] the volume of 
the fundamental cell is a significant fraction of the volume inside the surface of last 
scattering. Therefore, only the largest angular scales are modified but not the fine 
structures which determine the elongation. In Fig. |8] the mean ellipticity E and 
the elongation e computed from 50 torus realisations are compared with the elon- 
gations of the infinite flat ACDM model. As expected the mean values are almost 
identical and their small differences are conflned within the corresponding bands 
of cosmic variance. Thus, the anisotropy of the torus model is too weak to allow a 
discrimination between the models using the elongation measures. 

In addition to the above investigations we also analysed our simulated sky maps 
at other thresholds v. Quantitatively our statements are equivalent for hot spots 
and cold spots which is due to the isotropy properties of the random field. So we 
restrict ourselves in this paper to hot spots only. Nevertheless, if a single realisation 
of a measured sky map is studied, both should be combined into one statistic for 
a better significance. Other thresholds ly with i^ ^ ±1 possess a smaller number of 
spots and, therefore, lead to larger variances which is counterproductive in order 
to distinguish between cosmological models. Besides the mean values E(iy, rindn) 
and e{i>, n. An) we also studied higher moments, higher central moments and their 
statistical interpretations. Since it turns out that the mean value provides the most 
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Fig. 8. The ensemble average of the mean eUipticity E(i/, n^Qi^) with u = 1 (left) and of the 
mean elongation e(v, n, An) with u = 1 and An = 20 pixels (right) is presented. The elongations 
E and e of the infinite ACDM model are shown as a solid curve and its cosmic variance as a grey 
band calculated from 1000 realisations. The elongations of the torus universe with L = 3.86 L^f 
are shown as error bars which represent the cosmic variance based on 50 sky simulations. In both 
cases the window function of the W-band (channel 4) of WMAP is taken into account, but the 
noise is neglected. The sky maps have a resolution of N^^^^ = 512. 

robust measure we restrict our discussion to them here. The moving average apphed 
in e(h', n, An) leads to the best discrimination between cosmological models, but the 
additional parameter An has to be chosen adequately, since too small values result 
in small spot numbers (lack on statistics) and too large values result in a blurring 
caused by large spots. 



6. Summary 

In this paper the structures of CMB sky maps are studied with respect to the crucial 
question whether these structures betray some information about the underlying 
cosmological model. There are various quantities to analyse the structures of the 
niveau lines and the focus is put on the elongation which can be described by 
the eUipticity E, Eq. ([9|), or the elongation measure e defined in Eq. (fTTj) . The 
elongation depends on the parameter a, Eq. ()14p . which in turn depends on the 
multipole spectrum Ci and thus on the cosmology. The theoretical dependence of 
various elongation measures on a is given for the case that the statistical properties 
of the temperature field are those of two-dimensional homogeneous and isotropic 
Gaussian random fields. Our analysis leads to the conclusion that this dependence, 
however, is weak. Since the resolution of the maps has a superior influence on the 
niveau lines, the results obtained from sky maps depend on the pixelization and on 
the beam profile of the detector. Both deteriorating restrictions are analysed. It is 
found that the elongation cannot be analysed without reference to the pixelization 
and the beam profile. 

The dependence of the elongation on the cosmological parameters can be in- 
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vestigated for niveau lines classified either by their temperature or by the area 
which they encompass. The curvature of the universe is revealed by the elonga- 
tion as demonstrated in Fig. [6l where the elongation is considered as a function 
of the size of the spots. The best discrimination between cosmological models is 
obtained by applying the moving average to e(;/, n, An). A resolution of at least 
o'fwhm = 20 arcmin is required in order to distinguish between cosmological models 
which differ in the curvature by AJ7tot = 0.05. This excludes the analysis of the 
ILC map of the WMAP team, whose resolution is three times lower. The required 
resolution is achieved by the W-band sky map which has, however, much stronger 
noise. The analysis of simulated maps with the noise properties of the W-band map 
shows that it is also not suitable for this investigation. The Planck mission will 
provide maps with a significantly higher resolution as well as lower noise such that 
there is the hope that the results of that mission can be analysed with respect to 
the structure properties of the CMB. 

We also study the elongation properties of a multi-connected space form, i. e. 
a model of the universe having a non-trivial topology, because it is not statisti- 
cally isotropic and thus violates the assumptions which have to be satisfied for the 
analytical expressions to be valid. No distinction is found between the trivial topol- 
ogy, i. e. the concordance model, and the non-trivial topology with respect to the 
elongation. This is expected since the topology modifies the physics on the largest 
scales and has thus only a modest influence on the elongation. At least as long as 
the fundamental cell is so large that it constitutes a significant fraction within the 
surface of last scattering. Thus the best prospects for elongation measures can be 
found in the detection of curvature. 

Appendix A. The Gaussian Random Field and the CMB 

In this Appendix, we outline the derivation of the peak density and the ellipticity 
of the temperature fiuctuations ST (n) of the CMB under the assumption that the 
temperature fluctuations behave as a homogeneous and isotropic Gaussian random 
held 

on the 2-dimensional sphere. Here CTq (n) := {ST (n) ST (h)) is the variance of the 
temperature fluctuations. In an isotropic model ctq is independent of the direction 
n and is given in that case by 



{21 






Ci\Fi\' = C{0) . (A.2) 



The ensemble average is specified as (...), Fi accounts for the appearance of a beam 
profile and C(0) is the correlation function of the temperature fiuctuations in a 
homogeneous and isotropic universe at a separation angle i? = 0. To use a compact 
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notation in this Appendix, the subscript Fi at the temperature field STp^ is omitted. 
This index stands for the smoothing operation (|6]). The following derivations and 
notations are based on Ref. 1 

The temperature fluctuation 5T (fi) close to the direction fi' on the 2-sphere is 
given by the Taylor series 

5T (n) = 5T {n') + [V ,6T (n)] .^., {x' - x'') 

+ 1^^^^^^-Mn^(,^^,'^)(,._,0) + ... . (A.3) 



Vi is the covariant derivative on the 2-sphere with respect to the coordinate x^ . In 
the sequel we will choose x^ = 9 and x^ = 4>. The corresponding line element is given 
by ds^ — dO^ + sin^ OdcjP' . The associated covariant metric is 711 = 1, 722 = sin^ 9 
and "fij = otherwise. Therefore the Christofi^el symbols are Tfj = T^i — f^^, 
r22 = — sin 61 cos and T^,j = otherwise. 

The Taylor series (jA.3[) can be written in a compact way by introducing the 
variables 5T = 5T{n), -q^ = rj^in) := Vi5T{fi) and Cij = Qj (n) := V^Vj6T{n). 
Calculating the symmetric correlation matrix for Gaussian random fields with re- 
spect to the ensemble average in a homogeneous and isotropic model, one obtains 

2 2 

{6T6T) = 0-2, (STrj,) = 0, (77^77^) = -^7,^-, (STQj) = - yTij, 

2 2 

iViCjs) = 0, {(ijCrs) = -^ [lijlrs + lislrj + Ijslir] + -^lijlrs (A. 4) 

with 

al:^y}^^Ci\Fi\'l{l + l) (A.5) 



^-^ An 



and 



Y.^^^c^\Fi\'i^-mi + m + ^) ■ (A.' 



Appendix A.l. The peak density of the CMB 

Now the formulae for the number of extrcma per solid angle arc derived from the 
density of extrema ricxt('T-) = X^n ^ {n ^ fi-p) for an isotropic and homogeneous Gaus- 
sian random temperature field 5T (fi) on the 2-sphere. The density of extrema turns 
out to be given by 

next(n) = 5(^(n)) |det(C(n))| (A.7) 

by using x^—x^p sa rjj (n) (C"") (fip)- Here hp denotes the directions of the extrema. 
We choose ^ = f without loss of generality. We transform Qj onto its principal 
coordinate system by a rotation with an angle 9, thus obtaining the diagonal form 
— diag(Ai, A2), ordered by |Ai| > IA2I. At maxima or minima of the temperature field 
the eigenvalues Ai and A2 have both positive or both negative values, respectively. 
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At saddle points one of these eigenvalues is positive and the other is negative. In the 
following only maxima and minima are considered. The eigenvalues of the Hessian 
contain slightly different information as the eigenvalues of the inertia tensor of the 
area enclosed by the contour lines which are used to calculate the ellipticity E in 
Section [5l Now with the eigenvalues Ai and A2, the following variables are defined 



A1 + A2 
0-2 a 



Ai — A2 



Ai — A2 



2(Ai + A2) 2a2ax ' 



(A., 



in terms of which the Hessian Qj reads 



Cii 

C22 



(T2X a 



2 

<J2X a 



1 + 2 e cos(26') 
1 - 2 e cos(2^) 



C12 — —(T2xae sin(26') , 



(A.9) 
(A.IO) 
(A.ll) 



where a := \/l -\ — ^ and e is the elongation, Eq. ((Tl]) . It should be noted that 
in Ref. [I] e is termed ellipticity. In general these new variables are restricted to 
9 € [0,7r], X S (— oo,c») and e S [0,cxd). In case of extrema the interval of the 
elongation is confined to e £ [0, i). 

The transformation of the volume element in ('-space is 



dCiid(!^22dCi2 — lo^cc X dxdedO 



(A.12) 



The probability distribution for the variables v 



6T 



ff, X, e, and 9 is given by 



P{h', ff, X, e, 9) de dw dxd9 d fj — cxp 

t'2 



w 



exp 



- + 4e^a^ x^ 



X exp 



T]' 



2 dw dx d9 d fj 

[ax] e de ■ 



/27r V2^ 7^ T^f^i 



(A.13) 



where the result is simplified by using the variable w 



yr^ 



and introducing 



the abbreviation 7 := ^ — . 7 is determined by the power spectrum. Here w and 

x are independent and normalised (x^) = 1, (w^) = 1. 

As a result of the restriction to positive eigenvalues Ai > A2 > 0, i.e. restricting 



to maxima, one obtains x G [0, 00) and e G [0, i). Using this, 9 
and 



*2 2(7^ 



1 



det(C) = icr|x2a2(i_4e2) 



(A.14) 
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we obtain for the mean differential density of maxima 
Nn,^^{v,x,e,6) = / dw' (frf' dx' de'l d9' P{v' ,rf' ,x' ,e\0') 



X 5(77") S{x' - x) S{e' - e) Sid ' - e)^alx'^ a^l - A e'^) 



(A.15) 



2 


r h ir^) (r 


ix) exp 


w' 


^3 9*^ y 
X exp 


C [L 4C j [C 

1-7"^ 

- - + 4e^a^ x^ 

. V2 / . 


2 



The integration over the orientation angle yields (w = w{x, v)) 

2 



-/Vmax(;^,a;, e) 



7r2 61*2 0^:2 



e (1 — 4e ) (ax) exp 



X exp 
An integration over v in the last expression results in 



- ( - + 4e^a^ x^ 



(A.16) 



^max(a;, e)^ -^{-) (ax)'^ e (l - 46^) exp 



1 



- ( - + Ae^a^ I x^ 



(A.17) 



where we have substituted the variable v hy w. A further integration over x yields 

6 6(1-462)0* 



A^max(e) 



6l*27r (1 + 8^262) 



2„2^2 



(A.18) 



Here we have used the substitution y = (5 + Ae^o?^ x"^ . A final integration over e 
leads to the total number of maxima per solid angle 



A^max = / d6A^max(e) 
"'0 



87r6' 



1:2-1) Vl + 2a2 + 1 



2Fi(l,-;3;-2a2 



A■K9*^^/l + 2q;2 



Af, 



BfcE [(«^-i)vrT2^+i]V3 

Vl + 2q;2 



(A.19) 



where N^^^ :— . „^„ r- is the limiting case derived in Ref. ili We achieved the first 

max 4^g.2y'3 o i— i 

representation of A^ax by using the substitution y — e^ and the integral (2.2.6.15) 
in Ref. [5T]for the hypergeometric function 2^1(0, b] e; z). 

Integrating Eq. (JA.16P with respect to x and using Eq. (2.3.15.3) in Ref. [21] or 
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Eq. (2.3.15.7) in Ref. HU we obtain two representations of A^maxli^, e), 



Nma.^{v, e) 



(l-72)'a4e(l-4e2) 



^2/)*2 IQ^2^,2 



e^a'^{l--l'^) + iy 



X exp 



y.D^ 



" 4(l-72)(8e2a2(l_72)4.i) 
v/(l-7')(8e2a2(l-72) + i) I 



27re l-4e-' a"" 



7r2 6**2 ^8e2a2 (1 - 72) + 1 



94 



with q = —- 



and p 






exp 



Ap 



exp 



erfc 



2(1-7^) 
V2^ 



(A.20) 



^^_j — '' - 2(1-7=) 

and erfc(a;) the complementary error function. 
Integrating Eq. (jA.16|) over e, one gets 



Da{x) is the parabolic cylinder function 



N^^^{v,x) 



47r2 6'*2yri:^ 



exp 



w 

T 



exp 



X 



/(2;) 



(A.21) 



with f[x) := exp — (ax) — 1 + (ax) . An integration of Eq. (jA.21[) over i^ or a: 
yields 

-2- 
^max(a;) = ,„ ,3 ^,^ exp 



(27r)t6'* 






/(:^) 



and 



with 



A^max(i^) 



G(i^,7,a) := 71/(1-7^) 



exp 



(27r)5 61*2 



'2(1-7") 



exp 






G{v,j,a) 



(A.22) 



(A.23) 



V2t(1-7^) 
+ (a2(l~7')-l+7V) 



1 - Icrfc , ^'^ 

2 W2(l-7^), 



exp 



l+2a=(l-7=) 



^2a2(l-72) + l 



1 - icrfc ^" 

2 I V2(l-72)(l + 2a2(l-72)) 



respectively. We have used the integral (2.3.15.7) in Ref. [H] and the relation 
erfc(— x) — 2 — erfc(x) in order to compute iVmax('^)- 



July 19, 2010 0:43 WSPC/INSTRUCTION FILE paper 



24 R. Aurich et al. 



The mean differential density of maxima depending on the eccentricity e := 
VI-A2/A1 = 2[e/(l + 2e)]i/2 (e e [0, 1)) is given by 

N^Ae) - ^^'^'('~^') ^ (A.24) 



i^'^+e^)' 



and depending on the ellipticity E := y^hjh. = [(1 + 2e)/(l - 2e)]i/2 {E e [1, 00)) 

by 

j^ ^E)^ ^^-'EHE^-^) M25) 

^^^*^[E^l + 2a^) + 2E^l~2a^) + {l + 2a^)]i 

We obtain the corresponding expressions for the mean differential densities of 
the minima from those of the maxima by replacing v by —v and x by —x, e.g. 
-^min(j^) = -^max(— J^)- Using this wc get the combined distribution 

1 r ly'^^ 

(iV„iax + Nnun) {i^) = , ,3 . . exp -— (A.26) 





1 


" i/2- 

2 




(27r)i0*2^"^ 
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{a\l-l' 


)-l + 7^t^ 




L 







exp 



1+2q^(1-7^) 



In addition to the formulae given in Ref. [1} we have specified here analytical 
expressions for the densities N,^i,^{x,e), Naie,^{x), N,^i,^{e), A'^maxle) and Nmi,x{E). 
It should be pointed out that the formulae in Ref. [1] are obtained by the leading 
term of the Laurent series of our formulae at 0^2 = 1. As discussed in Section |4l the 
parameter a contains information about the underlying cosmology. 

Appendix A. 2. The ellipticity in the CMB 

In the following the formulae for the moments of the ellipticity E, the eccentricity 
£ and the elongation e at local maxima of an isotropic and homogeneous Gaussian 
random temperature field on a 2-sphere are derived. In the case of local minima 
the resulting formulae are also valid. The elongation which is considered in this 
Appendix results from a Taylor expansion at local extrema and can be computed 
from the eigenvalues of the Hessian (jA.8[) . For this reason the elongation e contains 
only information from local maxima or local minima. The same is valid for the 
ellipticity E and the eccentricity e, because they are related to the elongation by 
Eqs. (j9|) and (fTO| . respectively. This is a difference to the definition of the ellipticity 
in Section [5] where the eigenvalues of the inertia tensor of the area enclosed by the 
contour lines are used to calculate the ellipticity. 

Using the mean differential densities of the maxima with respect to the various 
arguments derived in [Appendix A.l[ one can define the distribution of a variable u 
subject to the constraint parameter z, 

P(^|,) ^ %^^ , (A.27) 
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e.g. P{e\v) = iVmax(i^, e)/A'iiiax(i^) with the variable e and the parameter v. 

The expectation values of the moments of the ellipticity with respect to the 
distribution 



p,. _ iVmax(e) _ 24 e(l - Ae^)a^^/lT2a^ 



^max (1 + 8a2e2)* [{q^ _ 1)^1 + 2^2 + l] 

16e(l-4e2) 



(A.28) 



(l + 8a2e2)^ 2-Fi(l,f;3;-2a2) 
are given by 

B{x, y) is the beta function. We get P{e) by using Eqs. (|A.18p and (|A.19p and the 
moments by applying the integral (2.2.6.15) in Ref. [21] From the expression of the 
moments, the mean value {n — 1) results in 



2a2+3-|./i±^ln 



-^ 



4[(a2_i)Vl + 2a2 + l] 
= 0.197 - 0.041 {a^ - 1) + 0{{a^ - if) 



(A.30) 



and the second moment [n = 2) in 

1 ~ 3 VI + 2a2 + 4(^1 + 2a2 _ i) / W 

(e2> = iLi^^ (A.31) 

4[(a2-l)VTT2^+l] 

= 0.049 - 0.018 (a^ - 1) + 0{{a^ - if) 

of the ellipticity at maxima. 

Using Eqs. (|A.24|) and (IA.19P we obtain the distribution 

P{e) = ^--(") 



96aVl + 2a2£3(i_g2) 



^A 32) 

[(a2 _ 1)7T~P2^+ 1] (4 - 4£2 + (1 + 2a2)£4)l 

from which the mean value 

(e) = f de P{e) e 
Jo 

= 0.715 - 0.059 (a^ - 1) + 0{{a^ - 1)^) (A.33) 
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and the second moment 



/ deP{e)e^ 
Jo 



_ 2 (a2 + 1 - VYT2^ 

~ (a2- 1)\/TT2^+1 

= 0.536 - 0.083 (a^ - 1) + OUa^ - 1)^) 

of the eccentricity can be calculated. 

Using Eqs. (|A.25p and (|A.19p one gets the distribution 



PiE) 



J ''max 



(A.34) 



(A.35) 



96aVl + 2a2^3(^2 _ ^^^ 



[(a2 _ l)Vl + 2a2 + l] (1 + 2a2 
which leads to the mean value 



-2(l-2a2)£;2 + (i + 2a2)^4p 



(i^> 



/OO 
dEP{E)E 

= -^ - V3i^(V2*) 



-^^.(72. 



-2 



10 



V33 
0((a' - 1)') 



2V3 



3)X(\/2i) 



VeVs 



- 2 E{^2i) 



[a^ - 1) 



= 1.648 - 0.217 {a^ - 1) + 0{{a^ - if) 
and the second moment of the ellipticity E 

/OO 
dEP{E)E^ 



(A.36) 



4^2 



[(a2 - l)Vl + 2a2 + 1] (1 + 2a2) 



3.131 



(A.37) 



0.980 (a^ - 1) + 0{{a^ - if) . 

Here K{k) is the complete elliptic integral of the first kind, E{k) the complete 
elliptic integral of the second kind and k the modulus. 

All these moments of the ellipticity E, the elongation e and the eccentricity e 
depend on a. On the other hand a is determined by the angular power spectrum 
of the underlying model, see Eq. p4p . For this reason also the moments of the 
ellipticity E, the elongation e and the eccentricity e depend on the angular power 
spectrum of the model. This dependence is considered in Section 14.21 
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